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ABSTRACT 



We present Spitzer Space Telescope and complementary ground-based in- 
""q^^I frared observations of the recurrent nova RS Ophiuchi, obtained over the period 
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64-111 days after the 2006 eruption. The Spitzer IRS data show a rich emis- 
sion hne spectrum superimposed on a free-free continuum. The presence of fine 
structure and coronal infrared hues lead us to deduce that there are at least two 
temperatures (1.5 x 10^ K and 9 x 10^ K) in the ejecta/wind environment, and 
that the electron density in the 'cooler' region is 2.2 x 10^ cm~^. The deter- 
mination of elemental abundances is not straightforward but on the assumption 
that the Ne and O fine structure lines arise in the same volume of the ejecta, the 
0/Ne ratio is >0.6 by number. 

Subject headings: stars: individual (RS Oph) — novae, cataclysmic variables — 
binaries: symbiotic — binaries: close — infrared: stars 



Introduction 



RS Oph is a recurrent nova (RN) that erupted in 1898, 1933, 1958, 1967, 1985, and 
possibly 1907 and 1945. The system consists of a semi-detached binary syste m comprised 
of a r oche-lobe-filling g iant (RG) and a massive (~ 1.2 M©) white dwarf (WD; IShore et al. 
19961 : iFekel et al.ll2000l ). As in classical novae, th e eruption follows a thermonuclear runaway 
on the surface of the WD (IStarrfield et al.lll988l ). In the case of the RS Oph class of RNe , 
however, the ejected material runs into, and shocks, a dense RG wind ( iBode fc Kahnlll985l ). 



The 1985 eruption of RS Oph was, for the first time, the su bject of a multi-wavelength 
observational campaign, f rom the radio to the x-ray (iBode 1119871 ). The most recent eruption, 
on 2006 February 12.83 (IHirosawal l2006l . we take this to define the origin of time: t = 0) , 



trig gered an even more intensive campaign. Infrared (IR) observations by iDas et al.l (120061 ) 
and 



Evans et a. 



(lO'Brien et al. 



Osborne et al 



(120071 . hereafter Paper I) show evidence for the shock, s een also at radio 
2o3and x-ray JBode et all bood : ISokoloski et"aD bood : iNess et al.l I2OO7I : 



20071 ) wavelengths, as the ejecta interact with the RG wind. HST observations 



of the ejecta from the 2006 eruption (IBode et al.ll2007l ) showed that emission in [Olll]5007A 
consists of two 'rings' consistent with coii stant velocity expansion of the radio lobes seen 
earlier in the eruption ( O'Brien et al.ll2006 ). In contrast the [Ne v]3426A lines possibly arise 
from two 'polar caps' at the extremities of the rings. 

We present he re further IR observations, obtaine d with the Infrared Spectrometer (IRS ; 
Houck et al.l |2004| ) on the Spitzer Space Telescope (I Werner et al.l l2003l : iGehrz et al.l 120071 ) 
and from ground-based facilities. 
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2. Observations 
2.1. Spitzer 

RS Oph was observed on 2006 Apr 16 {t = 62.51 days) and 2006 Apr 26 {t = 72.54 days) 
with SpitzerIRS as part of the Directors Discretionary Time program PID 270. Observations 
were performed using all IRS modules, with IRS blue (13.3-18.6 /im) peak-up on RS Oph. 
The spectroscopy consisted of 5 cycles of 14 second ramps in short-low mode, 5 cycles of 
30 second ramps in both short-high and long-low modes, and 5 cycles of 60 second ramps 
in long-high mode. IRS basic calibrated data products (BCDs) were processed with version 
14.0.0 of the IRS pipeline. Details of the calibration and raw data processing are specified 
in the IRS Pipeline Description Document^, version 1.0. 

Bad pixels were interpolated in individual BCDs using bad pixel masks provided by the 
Spitzer Science Center (SSC). For the low resolution data, multiple data collection events 
were obtained at two different positions on the slit using Spitzer 's nod function. The low 
resolution two-dimensional BCDs were differenced to remove the b ackground fl ux contri- 



bution and the data were extracted with SPICE (version 1.3-betal; ISpicd l2005l ) using the 
default point source extraction widths. For the high resolution data, the point spread func- 
tion nearly fllls the high resolution slit length, so it is not possible to perform background 
subtraction using nod pairs. Instead, separate sky observations were used to construct a 
master sky that was subtracted from the individual BCDs to remove the background con- 
tribution. Extraction was performed in the same manner as for the low resolution modules. 
For both resolution regimes the extracted, background- corrected data were combined, using 
a weighted linear mean, into a single output data flle. The point-to-point errors were es- 
timated from the standard deviation of the flux at each wavelength bin except when there 
were fewer than three data points in which case the errors generated by the SSC pipeline 
were added in quadrature to determine the flnal error. For the high resolution modules, in 
the order overlap regions, the long wavelength edge of the orders were much less reliable 
than the short wavelength edges, and so only the short wavelength regions of overlap were 
retained, i.e. the data from the lower order. As suggested by the SSC, data outside the 
ranges 5.2-14.5 fim for the SL module, 9.9-19.6 /im for the SH module, and 14.0-38.0 /im 
and 18.7-37.0 yum for the LL and LH modules, respectively, were discarded as unreliable. 
The spectral lines were flt using a least squares Gaussian routine that flt the line center, 
line amplitude, continuum amplitude and the slope of the continuum. Individual segments 
of the IRS spectrum were shifted upwards or downwards by < 10% to ensure that segments 



^ssc. spitzer. caltech.edu/irs/dh/PDD.pdf 
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adjoined smoothly. The spectrum for 2006 April 16 is shown in Fig. [U the spectrum for 
2006 April 26 is not substantially different. 



2.2. Ground-based observations 



The SpeX/IRTF data were obtained on 2006 May 1 (t = 77.61 days) and June 3 
{t = 110.72 days), UT using a 0.8" x 15" slit and a 10" N-S nod for background cancellation. 
No chopping was performed and extinction corrections were not necessary because of the 



proxi mity of the calibrator star. Data reduction was done using SpeXTools (jCushing et al. 



20M) with HD16471 6 (B9V) as the calibrator star. The flux of HD164716 was obtained 
by using the Kurucz (jl99ll : Il994j ) model of a Lyr scaled to the V magnitude of HD164716. 
SpeXTools makes an automatic, and in this case small, correction for extinction based on 
the measured {B — V) colors. UKIRT observations were performed on 2006 April 16 and 24 
{t = days 62.72 and 70.66 respectively) using the UIST instrument; details of the observations 
and data reduction are given in Paper I. The spectra from both telescopes, in the 1-2.6 /zm 
and 2.8-4.0 /im ranges, are shown in Figs. [2l 

Optical {BVr 'i'z' ) photometry was obtained by the robotic 2m Liverpool Telescope (LT; 
Steele et al.ll2004j ). sited on the island of La Palma, Canaries. Observations commenced as 
soon as RS Oph had declined to a level that would not saturate the detectors; standard 
photometric reduction techniques were used. 



3. Discussion 

The IR spectrum contains numerous emission lines. In the Spitzer wavelength range 
these lines are superimposed on a continuum that declines with increasing wavelength as fx oc 
A~^, consistent with free-free emission. The emission lines include hydrogen recombination 
lines, together with fine structure lines, many of which are coronal; the hydrogen lines and 
free-free emission will be discussed elsewhere. We use flux ratios for lines from the same 
atomic species and within the same wavelength band to estimate the electron temperature 
in the region in which they originate; the line pairs are listed in Table [H In Paper I the 
[Sivi]/[Six] flux ratio enabled us to estimate the temperature of the shocked, IR-emitting 
gas, to be ~ 9.3 x 10^ K. 

In this paper we have a wide range of ionic species and ionization stages at our disposal 
from which to estimate the electron temperature. Table [T] summarizes the collision strengths 
[Vt) at 10^ K (Paper I) for the species and ionization stages observed in the IR spectra. This 
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temperature is the highest in these sources available in common to all species and stages of 
ionizatio n, and the temperature-depe ndence of Q is relatively weak. The ionization fractions 
are from [Sutherland fc Dopital (119931 ) . The [Mg vil]5.50 fxm feature is blended with the H 16- 
7 recombination line (A = 5.52 /xm) at the resolution of the Spitzer SL mode. To correct 
for this we assume Case B (jOsterbrock fc Ferlandll2006l ). to estimate the expected flux ratio 
/(H16-7)//(H9-7) ~ 0.19 at electron density 2 x 10^ cm~^ (see below); we use the measured 
flux in the H 9-7 line (A = 11.31 /im) to correct the [Mgvil] line for the blend with H 16-7. 
In view of the uncertainties (particularly in the applicability of Case B) we should give the 
temperature derived from the Mg lines lower weight than that derived from other species. 

The derived temperatures are given in Tabled] (column 4). The uncertainties in logT(K) 
arising from uncertainties in the line fluxes are ±0.2. We note that the temperatures in 
Table [U derived from the Ne lines (mean ~ 1.5 x 10^ K), are lower than those derived from 
Si and Mg lines. The temperature derived from the Si lines in this paper is not significantly 
different from that given in Paper I during the early {t = 55.7 days) phase. Apparently there 
is a range of T values present in the shocked region. This result, derived from analysis of the 



IR em ission lines, is also suggested by x-ray data obtained on 200 6 April 16 (Osborne et al. 



20061 ) and follows naturally from models of the shock in RS Oph (jO'Brien et al.lll992l ). 



We shall assume that the electron temperature in the region in which the various Ne and 
[O iv] lines originate is 1.5 x 10^ K, while that in the Si/Mg line emitting region is 9 x 10^ K. 
The resulting critical densities ric, below which radiative de-excitation begins to dominate 
over coUisional de-excitation at the assumed temperatures, are given in Table [2J 

We use the fluxes in the [Ne v] 14.3,24.3 /xm lines to estimate the electron density in the 
'cooler' line-emiting region. Denoting the ^Pa, ^p^ and ^Pq levels as levels '2', '1' and '0' 
respectively, the fluxes / in the 14.32 /im and 24.3 /im lines are 

V 



f - A 



1 



where Uu (in cm 



u 



2, 1) is the population of the upper level, A is the wavelength, A 



(s~^) is the appropriate Einstein coefficient, V is the emitting volume and D (= 1.6 kpc; 
Bode 1119871 ) is the distance. Assuming that the [Nev] 14.2,24.3 /xm lines arise in the same 
region, the flux ratio gives ^2/^1 ~ 0.6. Detailed balance between radiative de-excitation, 
and electron collisional excitation and de-excitation amongst the three levels, together with 
the above values for n^jnx and the electron temperature, provide a value for the electron 
density n^. We find ~ 2.2 x 10^ cm~^ in the [Nev] emitting region. 



On the basis of radio imaging of the ejecta, lO'Brien et al.l (120061 ) estimate that the 
density in the radio-emitting shell was ~ 10~^^. K crn ~^ at t ^ 14 days. The wind/ejecta 



geometry is known to be complex (IBode et al.l 120071 ): however assuming that the density 



- 6 - 



declines as (appropriate for a uniformly expanding shell of constant thickness), we expect 
the electron density would be ~ 3 x 10^ cm~^ at t ~ 63 days, consistent with our estimate 
from the Ne lines. From Table [2] we see that this is generally below the critical density for 
all Ne lines with the exception of [Ne v]24.3/xm, the weakness of which compared with the 
other Ne lines may be consistent with a degree of coUisional de-excitation. 

We estimate the mass of emitting Ne, assuming that no + +?^2 = n([Ne v]), and using 
the [Nev] fluxes from Table [2] and equation We find n([Nev]) ~ 8.4 x IO'^^Mq for a 
distance of D = 1.6 kpc and, using the fractional abundance of [Nev] at 1.5 x 10^ K from 
Sutherland & Dopita (1993), we get M(Ne) ~ 2 x 10"^ {D/kpcf M©. In principle, we could 
also determine the mass of O; however it is clear from Table [2] that the deduced rie in the 
Ne/0 region exceeds the Uc for the ^P3/2 state of [O iv]. Consequently we conclude only that 
M(0) > 1 X 10-^ (D/kpc)2 Mq. 

In principle it would be straightforward to use the fluxes in the IR Ne (and other 
species) lines, together with the fluxes in the H lines, to estimate the Ne:H ratio in the 
ejecta /wind mix. However given the c omplex geornetry of the ejecta as evidenced by the 
radio (j O'Brien et al.l |2006| ) and HST (IBode et al.l 120071 ) observations we do not consider 
that this is justified until we have a better understanding of the environment of RS Oph. 
However assuming that the emitting volur nes of Ne and O coin cide, we find n{0)/n(Ne) ^ 0.6, 
compared with a {n(O)/ra(Ne)}0 of 6.6 (lAsplund et al.l 120051 ). 



In Fig. [3] we plot the time-dependence of the coronal line fluxes for which we have the 
longest time-base, namely S and Si; line fluxes for earlier data are taken from Paper I. It is 
curious that the S line fluxes decline monotonically over the first ~ 80 days of the eruption 
whereas there is a distinct minimum in the line fluxes around t ~ 70 days for the Si lines 
(the errors in the measured fluxes are typically ~ 5%). We consider that this result is real, 
because (a) there is a consistent pattern in the behavior of the Si and S lines and (b) it 
coincides (within the time resolution) with a distinct 'kink' in the BVr'i'z' light curves. 



in either 


the X- 


Ness et al. 


2007) 



(lOsborne et al.l 120071 . see Fig. [31 which shows only the V and B light curves for clarity). 



In classical novae, fluctuations in the visual light curve are often related to fluctuations 
in the mass-loss from the WD and the behavior of the pseudophotosphere; however, it is 
difficult to see how this can also lead to changes in the coronal line fluxes. Accordingly we 
consider it unlikely that this behavior is directly related to the mass-loss. 



For the 1985 eruption of RS Oph, lO'Brien et al.l (119921 ) estimated that breakout of the 
shock from the edge of the RG wind occurred ~ 60-70 days after the eruption, corresponding 
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to 70-80 days after the 2006 eruption, allowing for the greater inter-eruption interval. While 
this is close to the day 70 event in the 2006 eruption, it is unlikely that the behavior seen 
in Fig. [3] is related to any simple break out phenomeno i i, wh ich will be far more complex 
than the spherically symmetric model of lO'Brien et al.l (119921 ). Further, while the electron 
density we determine is for the Ne/0 region, it is likely that this value is not atypical and the 
disparity between the ric (Table [2]) and rie values suggest that the behavior depicted in Fig. [3] 
is not a colhsional de-excitation effect. Fig. [3] suggests we may be seeing the combined effects 
of recombination, abundance gradients and element segretation, but a detailed discussion is 
beyond the scope of this paper. 



4. Concluding remarks 

We have presented the Spitzer IRS spectrum of RS Oph during its 2006 eruption, the 
first mid-IR observations of a RN in outburst. The IR spectrum, from 1 /im to 30 /im, 
shows the rich fine structure (including) coronal and H recombination line spectrum of the 
shocked RG wind. There are at least two temperature regimes in the shocked wind, and the 
deduced electron density is consistent with extrapolation from radio observations earlier in 
the eruption. 
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Table 1: Estimated temperatures from coronal lines. 



Line ratio 


Flux ratio 


n at 10^ 


logT(K) 


[Si VI] 1. 96 /im/ [Six] 1.43 /im 


0.26 


0.43, 0.13 


5.96 


[Mg vil]5.50 /im/[Mg v]5.61 fxm 


3.25 


1.36, 1.06 


5.85 


[Ne lll]15.55 /im/[Ne v] 14.32 /im 


1.18 


0.79, 2.16 


5.18 


[Nelll]15.55 /im/[Nev]24.32/im 


4.53 


0.79, 0.75 


5.15 


[Ne Vl]7.65 /im/[Ne ill] 15.55 /xm 


7.29 


1.70, 0.79 


5.34 


[Ne II] 12.81 /xm/ [Ne v] 14.32 /im 


0.99 


0.40, 2.16 


5.08 



"From iBerrington. Saraph k Tullvl JlQQSh CfSivil). Sutler fc Zeippei] Jl994 ([Mgv], [ NeniD. I Griffin et al 



(|200l[ ) ([ NenD. iGriffin fc Badnelll (120001) ([ Nevn. iLennon fc Blaid (|l994f) ([Mgvn]), iMitnik et all (|2001f ) 
r[Nevi]l Ehang. Graziani fc Pradhanl (|l994h ([Six]). 



Table 2: Critical densities 



Identification 


Term 


A(/im) 


ric (cm ^) ^ 


[Six] 


2po 
^3/2 


2po 
^1/2 


1.4290 


1.0 X 10^° 


[Si VI] 


2po 
^1/2 


2po 
^3/2 


1.9625 


1.7 X 10^ 


[Si VII] 


^Po 


-^P2 


2.4822 


7.0 X 10® 


[Mgviii] 


2po 
^3/2 


2po 
^1/2 


3.0269 


2.6 X 10® 


[Six] 


^Po 


-^Pl 


3.7495 


4.5 X 10® 


[Neil] 


^Pl/2 


-'Pl/3 


12.808 


2.8 X 10^ 


[Nev] 


^P2 


-^Pl 


14.311 


1.4 X 10^ 


[Ne III] 


^Pl 


-^P2 


15.550 


8.1 X 10^ 


[Nev] 


^Pl 


-^Po 


24.300 


2.4 X 10^ 


[Oiv] 


^P3/2 


-^Pl/2 


25.871 


2.0 X 10^ 



°At 1.5 X 10^ K, for Ne/0 lines, at 9 x 10^ K for Si/Mg lines; see text. 
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Fig. 1. — Spitzcr spectrum of RS Oph for 2006 April 16 UT. Principal H recombination 
{n — m) and fine structure lines are identified; many of the lines without identification are 
also H recombination lines originating in high [n > 15) levels. 
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Fig. 2. — Top: Ground-based (UKIRT and IRTF) observations of RS Oph, obtained on 
2006 April 24 (UKIRT) and May 1 (IRTF); UKIRT spectrum has been displaced upwards 
by 0.1 X 10^^^ W m~^yum~^. Coronal lines are identified. Bottom: Ground-based (UKIRTl, 
UKIRT2 and IRTF) observations of RS Oph, obtained on 2006 April 16 (UKIRTl), April 24 
(UKIRT2) and May 1 (IRTF); spectra UKIRT2 and IRTF have been displaced downwards 
by 0.1 X 10~^^ W m~^//m~^ and 0.2 x 10~^^ W m~^//m~^ respectively. Coronal lines are 
identified. Many of the narrow lines without identification are H recombination lines. All 
dates UT. 
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Fig. 3. — Decline in coronal line fluxes (• = [Sivi], o = [Six], □ = [Svill], A = [Six], ■ = 
[Svill]); ordinate (left axis) is line flux in 10"-^^ W m~^. Swift fluxes in 0.3-10 keV channel 
(flUed points); ordinate (left axis) is in counts s~^. Liverpool Telescope (BV) light curves 
(open triangles; right axis); B data have been displaced upwards by 1.5 mag for clarity. 
Straight lines are least squares fits to light curves for the 40 < t (d) < 75, 75 < t (d) < 120 
ranges, illustrating break in the light curve around t ~ 70 days. See text for discussion. 




